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ABSTRACT 

In a ACDM cosmology, the Milky Way accretes satellites into the stellar disc. We use 
cosmological simulations to assess the frequency of near disc plane and higher inclina- 
tion accretion events, and collisionless simulations of satellite mergers to quantify the 
final state of the accreted material and the effect on the thin disc. 

On average, a Milky Way-sized galaxy has 3 subhalos with w max > 80km/s; 7 
with w max > 60km/s; and 15 with "u max > 40km/s merge at redshift z £ 1. Assuming 
isotropic accretion, a third of these merge at an impact angle 9 < 20° and are dragged 
into the disc plane by dynamical friction. Their accreted stars and dark matter settle 
into a thick disc. The stellar thick disc qualitatively reproduces the observed thick 
disc at the solar neighbourhood, but is less massive by a factor ^2 — 10. The dark 
matter disc contributes pddisc = 0.25 — 1 phalo at the solar position. Although not 
likely to be dynamically interesting, the dark disc has important implications for the 
direct detection of dark matter because of its low velocity with respect to the Earth. 

Higher inclination encounters 9 > 20° are twice as likely as low inclination ones. 
These lead to structures that close ly resemble the inner/outer stellar halos recently 
discovered by ICarollo et alj (|2007t ). They also do more damage to the Milky Way 
stellar disc creating a more pronounced flare, and warp; both long-lived and consistent 
with current observations. The most massive mergers (v max ~ 80km/s) heat the thin 
disc enough to produce a thick disc. These heated thin disc stars are essential for 
obtaining a thick disc as massive as that seen in the Milky Way; they likely comprise 
some ~ 50 — 90% of the thick disc stars. The Milky Way thin disc must reform from 
fresh gas after z = 1. 

Only one in four of our sample Milky Way halos experiences mergers massive and 
late enough to fully destroy the thin disc. We conclude that thick, thin and dark discs 
occur naturally within a ACDM cosmology. 
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1 INTRODUCTION 

The energy density of the Universe is primarily composed 
of a cosmologic al constant A (72%) and cold dark m atter 
(22%) (ACDM; ISpergel et al.ll2007l ; ISeliak et al.ll2006l ). Cos- 
mological N-body simulations make accurate predictions for 
the evolution of the dark matter component, finding self sim- 
ilar dark matter halos that merge hierarchically, resulting in 
a complex phase space distribution criss-cro ssed by streams 
l|Moore et al.lll99d ; IVogelsberger et al.ll2008l ). 

The Milky Way provides a natural laboratory for test- 
ing such predictions. The Milky Way is composed of a stel- 
lar disc, an old stella r halo, a central bar/bulge of stars, and 
a dark matter halo ( Fre eman fc Bland-Hawthorrj [20021 ) . as 
summarised in Table [T] Many of these components separate 
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into finer substructures that are likely evidence of past ac- 
cretions. The most visible of these are stellar streams: the 
Sagittarius stream that completes nearly two wraps in the 
plane perpendicular t o the Galactic disc (jlbata et al.lll994l ; 
iMaiewski et alj|2003l) ; the 'Monoce r os' ring that surr ounds 
the Galactic disc (|lbata et al.ll2003l ; IConn et al.ll2007l ); and 
many smaller streams that have rece ntly been found in the 
Galactic halo (|Belokurov et alJliooil ). 

More subtle accretion signatures may also be found in 
the Milky Way disc and stellar halo. ICarollo et"al1 (|2007l ) 
find an inner and outer stellar halo that are kinemati- 
cally distinct, suggesti ve of separate accretion events, while 
iGilmore fc Reid (119831 ) were the first to find that the Milky 
Way disc separates into thin and thick components. The 
thick disc is older (> 8Gyrs), hotter [a z ~ 42km/s), 
more metal poor ({[Fe/H]) ~ —0.6), and more slowly ro- 
tating (by ~ 40km/s) than the thin disc (see Table [TJ. 
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Its age separation from the thin disc corresponds to red- 
shift z = 1 in ACDM, which is the epoch when the 
merger rate is calculated and observed to drop rapidly as 
the Uni verse changes from matter t o dark energy dom- 
ination llDalcanton fe Bernstein! 20021; iBensbv et all 120071 ; 



iGottlober et al.l l200ll ; iKampczvk et alj|2006l ) This is cer 
tainly indicative of a merger origin for the thick disc, though 
there may not be any single well-defined formation epoch 
i|Nordstrom et all 12004 fl Finally, the Milky Way thin disc 
presents evide nce for past acc retions in its spiral arms, warp, 
flare and bar |Gerhardll2002l ) - all of which may be induced 
by satellite-disc interac tions (see e.g. iDubinski et al.l 120081 ; 
IKazantzidis et al . 1 120071 ) . 

The Milky Way is not alone in these properties. 
Streams are now being found in extra-galactic sys- 
tems - mos t nota b ly our nearest neighbour And romeda 
jlbata et alj 120071 ; iMartinez-Delgado et all 120081). while 
thick discs a p pear to be ubiquitous ( Tsikoudi 



tmck discs a p pear to be ubiquitous (isikoudi 
Bursteinl [l979; Ivan Dokkum et al.l Il994l; iMatthewsl 



1977 



2000 



Dalcanton fc Bernsteinll2002l;|Pohlen et alj2004llbata et al 



Dalcanton 



(2005 ); lYoachim fc Dalcantonl|2005l ; lYoachim 
120061 ; lElmegreen fc Elmegreenl (20061 ). 

In this paper, we study how the Milky Way disc 
affects the accretion of satellite galaxies in a ACDM 
cosmology, and how these satellites in turn affect the 
Milky Way disc. The Milky Way disc is the dominant 
mass component of the Milky Way interior to the solar 
circle. It is important because dynamical friction against 
the disc causes sa tellites to be preferentially dragged into 
the disc plane (|Quinn fc Goodman] 1 19861 ; [Q uinn et al.l 
1993). As satellites are torn apart by tidal forces, they 
depos it both the ir stars and their dark matter into a thick 
disc l|Lakel Il989h . The latter point is the key new idea 
presented in this work: a dark matter disc must form 
in a ACDM cosmology and we set out to quantify its 
mass and kinematic properties. At the same time, satellite 
accretions leave interesting morphological signatures in 
the stars and we quantify these to o. Low inclination 
events give an accret ed thick disc (jWalker et alj 1 19961 ; 
iPenarrubia et ahl 120061 ); high i nclination ones contribute 
to the Milky Way s tellar halo dBullock fc Johnston! [20051 ; 
iGauthier et al.l 120061); and both heat and dis t ort the Milky 
Way t hin disc (IQuinn fc Goodman! 1 19861; I Walker et al.l 
1996: Huang fc Carlberd 1 19971; ISellwood et ail 1 19981; 
Velazquez fc White! 1 19991; iFont et alj l200ll ; lAbadi et al 



20031 ; lArdi et all 120031; iBenson et all |2004|; IGauthier et a! 
20061; lHavashi fc Chib " 



20061 ; IKazantzidis et alj 120071: 



Villalobos fc Helmil 
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Satellite heating of the thin disc has been a sub- 
ject of particular int erest in many previous papers. 
IKazantzidis et ahl (|2007l ) were the first to be able to heat the 
thin disc enough to form a thick disc0. This is because they 
measured and used the high-redshift satellite distribution 
that is both more massive and more destructive (because 
the orbi ts are more radial) t han the surviving distribution 
(see also lDiemand et al.l[20 04) . Using the redshift zero satel- 



1 A connection with merger s for the forma tion of the thick disc is 
compelling, but not certain. iKroupal d2002h . for example, suggest 
an alternative mechanism where the thick disc forms from massive 
star clusters produced in a phase of very high star formation rate. 



lite d istribution will not give enough heating to form a thick 
disc (|Font et alj l200ll ; IDubinski et all [2008). Here we also 
calculate and u se the high-redshift satell ite distribution. We 
find similarly to lKazantzidis et alj (|2007l ) that only the most 
massive of these satellites can heat the thin disc enough to 
form a thick disc. In addition, we consider what happens to 
the material accreted in such encounters - both the stars 
and the dark matter, and we quantify the fraction of the 
thick disc that can be accreted. This is particularly relevant 
for extra-galactic systems, where one counter-rotating thick 
disc has been found that cannot have forme d from thin disc 
heating alone (|Yoachim fc Da lcanton 2005). 

Our strategy is to use a cosmological structure forma- 
tion simulation to determine the frequency of satellite-disc 
encounters, and collisionless simulations of satellite mergers 
to quantify the final state of both the accreted material and 
the stellar disc. We do not model the gaseous component 
of the disc, since we are primarily interested in collisionless 
mergers. Models for the Milky Way suggest that stars have 
continued to form at a constant rate since z = 1, yet the sur- 
face density of ga s likely remained constant with time due to 
fresh gas inflow dRocha-Pinto et alj|2000l ; iNaab fc Ostrikerl 
20061; I Just fc~J ahreiss 2003). The gas is important because it 
can cool and re-form a thin disc in the 8 Gyrs between 2=1 
and the present. But it is also important because, unlike 
the stellar disc that is continually heated by mergers, the 
gas disc cools and maintains a low vertical dispersion. This 
keeps the density in the disc plane high, producing more 
significant dynamical friction. As a result, our gas free sim- 
ulations likely under-estimate the effect of disc plane drag- 
ging. To compensate for this, we consider a range of disc 
masses from half of the present Milky Way disc mass up to 
its present mass. We will study the effect of gaseous Milky 
Way discs in future work. 

A complementary approach is to run full gas- 
hydrodynamical simulations of the formation of the Milky 
Way that include cooling, star formation and feedback 
physics. Such simulations remain a significant challenge 
both because of the computational cost of including gas 



is most important (see e.g. Thacker & Couchmanl 


200C ; 


Ricotti et alj|2002; 


Abadi et al. 2003; Robertson et al.l 


2004; 


Ceverino & Klvoin 


2007: Governato et alj 2007). DesDite 


these difficulties, 


Abadi et al.1 (20031) found a thick disc 



in their simulated Milky Way that formed from a mix 
of accreted a nd he a ted stars, wh i le in a seri es of papers, 
iBrook et al.1 l|2004l) . iBrook et al] (|2005l ) and iBrook et al.1 
1 20071 ) find thick discs form directly from heated gas in their 
simulations. Our approach is complementary to these stud- 
ies. Predictions for the evolution of the dark component rely 
only on getting the gravitation al clustering right and are now 
robust (|Heitmann et al.|l2007l ). We combine these with our 
knowledge of what galaxies look like in the Universe today to 
largely side-step the unknown physics of galaxy formation. 



2 In a study concurrent with our own, Villalobos & Hclmi (2008) 
also consider the formation of thick discs from satellite heating 
and manage to successfully form heated thick discs. Their focus 
is on finding unique morphological and kinematic signatures of 
such accretions; ours is on calculating the expected number and 
mass of mergers in ACDM and the effect these have on an average 
Milky Way galaxy. In this sense our studies are complementary. 
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This paper is organised as follows. In Sj2]we use a concor- 
dance cosmology simulation to assess the frequency of high 
and low inclination satellite accretions for a typical Milky 
Way galaxy. In fj3] we present a suite of collisionless simu- 
lations of disc-satellite mergers, where we vary the satellite 
impact angle, orbit and mass. Finally, in Sj3]we present our 
conclusions. 



2 COSMOLOGICAL SIMULATIONS: THE 
LIKELIHOOD OF DISC ENCOUNTERS 

In this section we use a cosmological simulation to assess 
the likelihood of near-disc plane mergers in the Milky Way. 
Some movies of the simulation and the subhalo orbits can 
be seen at: http://justinread.net 



2.1 Description of the simulation 

We use t he cosmological ACDM simulation already pre- 
sented in iDiemancT et al. (2005|). The simulation was run 
using PkdGRAV (Stadcl 2001), with cosmological param- 
eters: (fi m ,OA,<T8,ft)=(0.268,0.732,0.7,0.71), and a box of 
size Lbox = 90 Mpc, with 300 3 particle s. The initial cond i- 
tions were generated with GRAFIC2 i|Bertschingerl l200ll ). 
From the simulation volume, we extracted four Milky Way 
sized halos at a mass resolution of m p = 5.7 x 1O 5 M0. 

The subhalos inside each 'Milky Way' and at each 
redshift output were identified using the Affifl algorithm 
jGill et all 12001 ). First, all particles are grouped together 
within isodensity regions. Second, unbound particles are it- 
eratively removed until a bound structure is converged upon. 
We then assign particles to the smallest structure they ap- 
pear in so that each particle is counted only once. We con- 
sidered all (sub)halos with > 50 particles. 

Once a subhalo was identified, we found its position 
by centri ng on it with the ' shrinking sphere' method de- 
scribed in lRead et alj (1200 6a) : we used the peak of the ro- 
tation curv^3 of all bound particles about its mass centre as 
a proxy for mass. Each subhalo was then traced backwards 
in time, its progenitor being the subhalo at the previous 
redshift that contains the majority of its particles. To avoid 
ambiguities in this halo tracking, we ordered subhalos by 
mass so that subhalos were linked to their most massive 
progenitor not already claimed by a larger subhalo. A final 
complication can occur if two subhalos are about to merge. 
In this situation AHF sometimes over-estimates the mass 
of the smaller of the two. We dealt with this problem by 
searching for sudden spikes in mass at pericentre and re- 
moving these by assigning instead the mass found at the 
previous output time. 

We are interested in finding subhalos that merge with 
the disc at a given redshift. However, standard merger trees 
will not suffice since these define a merger as being when a 
subhalo enters the main halo. By this definition the LMC 
has already merged with the Milky Way! Instead, we de- 
fine a subhalo as merged if it: (i) has less than a tenth of 
its peak circular speed considered over all times; and (ii) 



has passed within r n 



of the main halo after z = 4.35. 



Our results are not sensitive to changes in the peak circu- 
lar speed fraction or the redshift cut, but are sensitive to 
Emerge, particularly for subhalos with peak circular speeds 
Umax ,$40 km/s. For this reason, we present results for a 



range of r„ 



25,50, 75 kpc. 



2.2 The likelihood of disc plane mergers 

Figure[T] left panel, shows the cumulative maximum circular 
speed function of surviving (dashed) and accreted (solid) 
subhalos for the four Milky Way-sized halos. There is some 
cosmic variance amongst the halos - in particular halo 1 is 
distinctly different from the other three. It has three very 
massive mergers with f ma x > 90km/s that occur at z < 1. 
As we will show in section 33 these are very damaging to a 
disc, while their late accretion leaves little time to rebuild a 
thin disc. As such, it seems likely that halos 2, 3 and 4 are 
more representative of the Milky Way. 

Our results f or the surviving subhalos agree well with 
previous studies (|Gao et alj 12004 ; iLudlow et al.l 120081 ) and 
have cumulative i> max functions well fit by their fitting func- 
tion: 



N(> w max ) 



host J 
' max / 



(i) 



where v^i* is the peak circular speed of the host galaxy; 
q = 2.559; and A = 0.06. (Note that this underestimates 
slightly the low mass end. We find a = 3 g i ves a better 
formal fit in agreement with iDiemand et al.l (|2004l ).) Our 
accreted subhalos are well fit by a similar function with an 
exponential cut-off: 



N(> v max ) = A 



\ 7 ,host 
^ ^max 



exp 



/ ^max \ ' 



(2) 



3 http : //www . aip . de/People/AKnebe/AMIGA/ 

4 Also called the 'circular speed curve'. 



The fitted values depend on our choice for the subhalos' 
point of closest approach, r mC rg C . Figure [1] shows results for 
Emerge = 50 kpc, for which we obtain best fit values: A = 1.4, 
vq — 0.63, a — 2, and (5 — 6. However, decreasing r mcrgc 
decreases the number of lower mass subhalo mergers. For 
'"merge = 25 kpc we find: A — 4.5, vo = 0.48, a = 1.36, and 
P — 2.77, while r mcrgo = 75 kpc gives: A = 1.3, v = 0.64, 
a = 2.1, and (3 = 6.6. We likely over-estimate the pericen- 
tres because of our limited number of simulation outputs, 
and because, if we were to include stars and gas in these 
simulati ons, the subhalos w ould be more resilient to tides 
(see e.g. iMaccio et al"1l2006T l. However, our results are well 
converged for w ma x > 40km/s and r- mcr gc < 75 kpc. Since it 
is these most massive halos that are of most interest for ei- 
ther damaging the disc or depositing accreted material, our 
choice of r m0 r g c is not critical. We assume r mer gc = 50 kpc 
from here on. 

From our above fits, assuming ifex = 200 km/s, we 
find: 3 mergers with w ma x > 80 km/s; 7 with w ma x > 60 km/s; 
and 15 with v ma x > 40 km/s. Nearly all mergers at t> m ax > 
60 km/s are complete by, or shortly after, z ~ 1 (see the mid- 
dle panel of Figure [TJ) . Most of the massive accreted subha- 
los were found to be on highly eccentric orbits, with e > 0.8 
(see right-most p a nel of Figu re [p. This agre e s wel l with 
iKazantzidis et al.l ([2007) and IDiemand et all |2004) who 
find that the accreted subhalos show a more eccentric orbit 
distribution than the survivors. It is important to stress that 
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Table 1. The distinct components of the Milky Way. From left to right the columns give: mass; local volume density; local surface 
density for z < 1.1 kpc; half mass scale length and height (see i|3.2ll ; R,<f>,z velocity dispersion; rotation velocity; mean metallicity; 
mean age; and references. Where two values are given this reflects the fact that parameters often depend on the magnitude cut, or in 
which direction one looks (see e.g. iKaraali et al.|[2004t ljuric et alJIgOQsl: ICabrera-Lavers et al1l2007h. Similar results have been observed 



f or nearby extra-galactic sy ste ms (see tex t). Data are taken from: l lKuiiken &: Gilmoreill989l; 2 Morrisorjlll993l; 3lDehnen fe Binnevl 
4 lWilkinson fc Evans! Il999l: 5lOihal l200ll: 6lSoubiran et al] 120031: 71 Karaali et al .1120041: 8lHolmberg fc Flvnnl l2004i ~9l Juric et all 
lO lFulbrieht et al.ll2006l ; ll lBensbv et alj|2007l ; 12 lCarollo et al.ll2007l; 13 lSeabroke fc Gilmorell2007l . 
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Figure 1. Left: cumulative maximum circular speed « ma x functions for accreted (solid lines) and surviving (dashed lines) subhalos. We 
include only halos with greater than 50 particles. For the accreted halos, we consider only halos that at z = have 0.1 of their peak 
» msx considered over all times, and which pass within r mer ge = 50 kpc of the disc after 2 = 4.35. Four host halos are shown, taken from 
the concordance cosmology simulation described in J2] (1 black, 2 red, 3 green, 4 blue); their v^ fj ar e marked in km/s. T he thick solid 
line shows a power-law fit to the surviving subhalo t) max function taken from lGao et al. I 1120041) and iLudlow et al. I [12008ft . A fit to the 
accreted subhalos is also shown (thick dotted line; see text for details). The vertical dotted lines mark u max = 40, 60, 80 km/s assuming 
u ma5 = 200 km/s. Middle & Right: evolutionary tracks in u max and radius r for the six most massive subhalos with u ma x > 40 km/s 
that are accreted by z = (solid, dotted, short dashed, dot dashed, triple-dot-dashed and long dashed lines). The different colours 
correspond to the same four host halos as in the left panel. For halo 1 (black lines) the u max are marked on the right panel in km/s. 
Redshift z = 1 (~ 8 Gyrs) is marked by the vertical solid lines. Note that the pericentres are likely over-estimated due to the limited 
number of simulation outputs. 



the accreted subhalos are systematically more massive than 
the su rvivors. As pointed out recently bv iKazantzidis et al.l 
(2007), it is essential to consider the accreted satellite dis- 
tribution when studying satellite-disc interacti ons, other- 
wise the effects will be underestimated (see also lSales et al.l 
2007). Here we add the additional point that it is impor- 
tant also to consider only the accreted halos that actually 
interact with the disc. Similar results have been derived pre- 
viously from semi-analytic models ( Taylor fc Babull 200 ll . 
12004 120051 ; IZentner fc Bullockll2003l ; Izentner et al]|2005h . 

Assuming isotropy, the probability that mergers will 
lie within an angle 6 of the disc is given by P = sin#. 
Hence, roughly a third of the mergers will occur with 
9 < 19.5°, a third with 19.5° < 6 < 41.5°, and a third 
with 41.5° < 8 < 90°. Several recent studies - that focus 
on the surviving satellites - have found anisotropic satel- 



lite distribu tions, but some find polar alignmen t (the Holm- 
berg effect; lHolmberdll974l;lBailin et al.ll2007h. others find 



planar alignment (the anti-Holmberg effe ct; iBrainerd 2005; 
ISales et al.ll20Q7| ; |Faltenbacher et aill2007h. and some find ! 



statistically significant alignment at all l Azzaro et al.|[2006l ). 
From our four sample halos, we find two that have a very 
anisotropic merger history; indeed infall along filaments is 
expected to give rise to some anisotropy dKnebe et al.ll2004l ; 
IZentner et al.1 2005; Libe skind et al.|[2005l ). However, we can- 
not address the issue directly with our simulations, since we 
do not know how the disc shoul d align with the halo (see 
also lAgustsson fc Brainerdl 120061 ) . In this work we will as- 
sume isotropy. 

Putting all of the above together, we find that a typical 
Milky Way sized halo will have 1 subhalo merge near the disc 
plane (0 < 20°) with f ma x > 80 km/s; 2 — 3 with i> ma x > 60; 
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and 5 with w ma x > 40 km/s. Away from the disc plane there 
will be twice as may mergers at the same mass. 

Our results a gree well with a recent study by 
IStewart et ail (|2007h . They use cosmological simulations to 
study a much larger sample of Milky Way halos than we do 
(17,000), but at lower resolution. They conclude that 95% 
of all Milky Way halos have more than one merger with a 
subhalo m > 5 x 10 10 M Q (u max ~ 75 km/s), while 70% have 
a merger with a subhalo m > 10 11 Mq (w max ~ 90 km/s). 
(They define a merger as being anything that enters the 
virial radius, which corresponds to the sum of our dashed 
and solid lines in Figure [T] left panel.) They suggest that 
these latter mergers could prove problematic for thin disc 
survival. In i|3] we will show that such mergers do signifi- 
cantly heat the thin disc, but they do not destroy it; indeed, 
these mergers (and those at slightly lower mass) are likely to 
be essential for forming thick discs as massive as that seen 
in the Milky Way. 

As a final point, we remark that it is surprising that 
subhalos with i; max ~ 60 — 70 km/s manage to m erge in 
significantly less than a Hubble time. [Zhao ( 2004|) showed 
recently that, once mass loss due to tides is included, dy- 
namical friction times for such subhalos in the Milky Way 
are greater than 1 Gyrs. Indeed, we wi ll verify that this is 
the case in j5X2] l|Gauthier et alj (2006) find similar results 
in their simulations). Yet such subhalos clearly do merge 

- and rapidly - in the cosmological simulations. The solu- 
tion is presented in Figure [2] where we show logarithmic 
density contours of a spherical region of radius 100 kpc se- 
lected around the most massive subhalo in halo 2 at z = 2.4. 
These same particles are tagged and plotted at later times, 
z = 1.9, z — 1.6, z = 1.4, and z = 1.2. Notice that, although 
this subhalo has Umax = 75 km/s (m = 2.55 x 10 10 Mq), a 
much larger region around the subhalo enclosing 10 1 Mq 
(and some other smaller subhalos) largely co-moves with it 

- these subhalos fall in inside a larger loosely bound group 
environment and are effectively deposited at lo w apocentre. 
Simila r results have been presented recently bv lLi fc Helmil 
l|2007i ). who suggest that this might explain the apparent 
alignment of many Local Group gal axie s in planes on t he 
sky (see also lLake fc D'Onghiall200sl and lD'Onghiall2008r ). 

To approximate the effect of the larger group environ- 
ment, for most of the simulations that follow in S|3] we will 
use a small initial apocentre of r a = 30 kpc. We justify this 
further in M3.3.2I where we show the effect of a larger apoc- 
entre orbhfj. 

The above is particularly interesting in the context of 
the Sagittar ius dwarf wh ich is currently merging with the 
Milky Way. IZhaol (|2004l ) pointed out that it would have 
to have been uncomfortably massive in the past in order 
to arrive at its present location. This tension is resolved 
if Sagittarius accreted within a larger group environment. 
Modelling the effect of this in detail could be essential for 
a proper understanding of the Sagittarius stream and the 
orbits of the Large and Small Magellanic Cloud^f] 

5 iHuane fc Carlberd <1997ft argue, quite reasonably, for starting 
satellites at larger apocentre since this should be more cosmologi- 
cally consistent. Here we suggest that it is actually less so because 
o f this larger group en vironment . 

6 iBerrier et al.l J2008f) have recently studied isolated versus group 
mergers in clusters and find that 70% of cluster galaxies are ac- 



3 CONTROLLED SIMULATIONS: DRY 
SATELLITE MERGERS 

In iJ21 we demonstrated that Milky Way sized galaxies will 
have on average 2 — 3 satellites merge with u max > 60 km/s 
near the disc plane (9 < 20°) at early times, and twice 
as many at higher inclinations. In this section, we use 
a suite of collisionless simulations to measure the mor- 
phology and kinematics of the material accreted in such 
mergers, and the effect on the Milky Way disc. Some 
movies of our reference simulation LMC-10°can be seen at: 
: // j ust inread . net 



3.1 Description of the simulations 

We set up our Milky Way model (disc+halo system) by adi- 
abatically growing a disc inside a spherical halo. The initial 
spherical halos were generate d from a distribution fun ction 
using the method described in Kazantz idis et al.l (|2004| ) with 
the added refinement that the halo is composed of two mass 
species arranged on shells. The outer shell had particles 5x 
more massive particles than the inner shell, in order to in- 
crease the effective resolution in the central parts. We then 
inserted a massless disc of particles with scale-length, Rd — 3 
kpc. We slowly grew the mass of this disc over 1.5 Gyr to a 
final mass of 3 x 10 10 M (MW) and 6 x 10 10 M (MWB) 
while holding the disc particles fixed in place. After this 
time, we set the kinematics of the disc to give a constant 
Q = 1.5. Evolved in isolation, the MW system formed a 
small, weak bar after a Hubble time but otherwise did not 
change significantly; by contrast, the more massive MWB 
formed a strong bar. Recall that we do not know the to- 
tal baryonic mass or gas fraction of the Milky Way disc at 
z = 1 when most of these mergers took place, so it is use- 
ful to consider a range of disc masses. Our MW model has 
a peak circular velocity of 220 km/s and the dark matter 
halo provides 200 km/s of this (see Figure [3}. Current mass 
models for t he Milky Way favou r lower halo contributions 
(~ 170 km/s; Iklypin et aj]|2002ft . so our MW model corre- 
sponds to a minimum disc situation; MWB examines the 
effect of a more dominant disc. 

We chose three models for our satellite: Fornax, LMC 
and LLMC. Fornax, with dark matter contribution to the 
peak circular velocity u max = 25 km/s, models the satel- 
lite on the Fornax dwarf spheroidal galaxy which has a 
visible mass of ~ 5 x 10 7 Mq, and a dynamical mass of 
M{< 1.5kpc) ~ 7.5 x 10 8 M Q (| Walker et alj l2006h . It was 
initialised as a two component spherical galaxy, with star 
and dark matter velocities dra wn from numerically calcu- 
lated distribution functions as in lRead et al.l (|2006al ) . LMC, 
with u m ax = 60 km/s, models the satellite on the Large 
Magellanic Cloud - the largest surviving satellite of the 
Milky Way - with visible mass of M ~ 2.7 x 10 9 M Q , and 
a dynamical mass of M(< 8.9 kpc) = 8.7 ± 4.3 x 10 9 M Q 
l|van der Marel et alJl2002T ). LLMC, with u max = 80 km/s, is 

creted in isolation. However, their study focuses on the surviving 
subhalos rather than the accreted ones, and on cluster environ- 
ments, rather than the galaxy environment we study here, so 
detailed comparisons are necessarily difficult. We will pursue this 
issue of group versus isolated mergers in more detail in future 
work. 



http 
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Figure 2. Subhalos fall in inside loosely bound larger groups. Shown here are logarithmic density contours of a spherical region of radius 
100 kpc selected around the most massive subhalo in halo 2 at z = 2.4, before this subhalo falls into the host halo. These same particles 
are tagged and plotted at later times, z = 1.9, z = 1.6, z = 1.4, and z = 1.2. The peak circular velocity of this subhalo is 75km/s, while 
its spherical overdensity mass and radius at z = 2.4 are Mgo = 2.55 X 1O 1O M0 and rgo = 15.8 kpc, respectively. The much larger region 
extending to 100 kpc and enclosing 10 11 Mq (and some other smaller subhalos) largely co-moves with this subhalo. 



a 'large' LMC model that considers the effect of the most 
massive mergers found in Sj2] Both LMC and LLMC were 
set up as scaled versions of our MW Milky Way model. Our 
model satellite properties are given in Table [2] and Figure [3] 
The total peak circular speeds u^ax f° r Fornax, LMC and 
LLMC are t)^ = 30, 70, 90km/s, respectively. These are 
all larger than the quoted « max above because those were 
for only the dark matter contribution to the peak circular 
speed. We quote v max rather than D max from here on for 
easier comparison with the results from ij2] 

For the orbit parameters, we chose a wide range of ini- 
tial inclination angles to the disc from 10 — 60°, one retro- 
grade orbit (Ret-10°), and range of pericentres and apocen- 
tres, as detailed in Table [2] We chose a typical apocentre 
of 30 kpc for reasons outlined in 33 We test the effect of a 
larger apocentre in one run: LMC-10°-abig. 

The simulations were evolved using the collisionless 
tree-code, PkdGRAV (|Stadel|[200ll '). We used a maximum 
timestep of 0.005 Gyrs with nine rungs of decreasing size 
in powers of two. The final evolved systems were mass and 
momentu m centred us i ng the 'shrinking sphere' method de- 
scribed in Re ad et"al1 (|2006al ). and rotated into their mo- 
ment of inertia eigenframe with the z axis perpendicular to 
the disc. 



3.2 Some definitions and nomenclature 

Some of the density profiles we obtain are not easily fit by 
a simple functional form. For this reason we use half mass 
scale length^ as this is then model independent. Discs are 
typically fit with an exponential profile, given by: 



p(R,z) 



M 



4ttR%z 



-R/R 0e -\z\/z 



(3) 



7 There can be some confusion here as 'half mass scale length' 
depends on the choice of coordinates. We work in cylindrical polar 
coordinates (R,ij>,z) so that R\/2 means the radius containing 
half of the mass, having already integrated over z and <j>- Similarly 
Z1/2 means the \z\ that contains half of the mass, having already 
integrated over R and <j>. In spherical polar coordinates (r,8,rfr), 
the half mass radius rj/j is the radius containing half of the mass, 
having already integrated over 9 and <f>. So rj/ 2 Ri/2 even in 
the limit of spherical symmetry. 
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Figure 3. Rotation curves for the different galaxy mod- 
els. The upper red line is for MWB. The dotted lines show 
the dark matter halo contribution . The data point s show 
the mean of HI measurements from iBania fc Lock manl J19841). 
IWeaver fc Williams! lTl974h,IWeaver fe Williams! lll973l) , lMaIhotral 
Jl995h and lKerr et al.l lll986h . The blue dashed line shows the ro- 
tation curve from one of the subhalos taken from the cosmological 
simulations of ij2] 



where Mo is the disc mass, and Ro,zo are the scale length 
and height, respectively. For this distribution, the half mass 
scale lengths are related to the scale lengths as: R1/2 = 
1.68-Ro, Z1/2 = 0.7zo • Other common variants have p(z) oc 
sech 2 (z/2 s ) and p(z) oc e~ < - z / z "' 1 , whose scale lengths inter- 
relate as: Z1/2 — 0.55z s = 0.48z e = 0.7zo- We will refer only 
to the half mass scale lengths from here on. 

We also refer frequently to the 'solar neighbourhood'. 
By this we mean a cylindrical slice centred on Rq — 8.5 kpc, 
lkpc thick - i.e. 8 < R < 9 kpc; a slice 7.5 < R < 8.5 kpc 
gives nearly indistinguishable results. 
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Galaxy model 


Description 


N* 




e* (kpc) 


eDM(kpc) 


M»(M ) 


r*(kpc) 


^dm(M ) r DM (kpc) 


MW 


Milky Way; light disc 


7.5 x 10 5 


2 x 10 6 


0.06 


0.1; 1 


3 x 10 10 


3 


10 12 25 


MWB 


Milky Way 


7.5 x 10 5 


2 x 10 6 


0.06 


0.1; 1 


6 x 10 10 


3 


10 12 25 


LLMC 


Large-LMC 


7.5 x 10 5 


2 x 10 6 


0.015 


0.024; 0.234 


3 x 10 9 


1.78 


1 x 10 11 15 


LMC 


LMC-like 


7.5 x 10 5 


2 x 10 6 


0.015 


0.024; 0.234 


7 x 10 s 


0.7 


2.4 x 10 10 6 


Fornax 


Fornax-like 


10 e 


10 6 


0.005 


0.02 


5.6 x 10 7 


0.1 


10 9 2 


OVYYlUlCL~t%OTl 


(x, y, z) (kpc) 


{Vx ,Vy,v z ) (km/ s ) 






notation 


-^outv^yrsj 






LMC-10 


(29.5,0.27,-5.2) 


(-6.3,89.3,0.35) 


10° 


0.8 


prograde 


6.5 






LMC-20 


(28.2,0.12,-10.3) 


(-2.2,80.1,0.82) 


20° 


0.8 


prograde 


4.3 






LMC-40° 


(23,0.12,-19.3) 


(-1.8,80.1,1.52) 


40° 


0.8 


prograde 


6.3 






LMC-6O 


(15,0.12,-26) 


(-1.2,80.1,2.0) 


60° 


0.8 


prograde 


7 






LMC-Rct-10° 


(29.5,0.237,-5.2) 


(-6.3,-89.3,0.35) 


10° 


0.8 


retrograde 


5.3 






LMC-10°-eless 


(29.5,0.27,-5.2) 


(-6.3,143,0.35) 


10° 


0.36 


prograde 


1.1 






LMC-10°-abig 


(80,0.27,-15.2) 


(-6.3,62.5,0.35) 


10° 


0.74 


prograde 


6.8 






Fornax-10° 


(29.5,0.27,-5.2) 


(-6.3,89.3,0.35) 


10° 


0.8 


prograde 


5 






LLMC- 10° 


(29.5,0.27,-5.2) 


(-6.3,89.3,0.35) 


10° 


0.8 


prograde 


5.5 






LLMC-10°-MWB 


(29.5,0.27,-5.2) 


(-6.3,89.3,0.35) 


10° 


0.8 


prograde 


4.5 







Table 2. Simulation labels and parameters. The top table shows the different galaxy models we use, two for the host galaxy (MW/MWB), 
and three different satellite galaxy models: LMC, LLMC and Fornax. The columns show from left to right: the number of star and dark 
matter particles, A^, TVdm; the force softenings, £*,£dm; and the masses and scale lengths of each component, M,,r,,MDM,DM. Where 
two force softenings are given, these are for the low mass inner and higher mass outer halo particles (see ^13. II for details). The bottom 
table shows the different initial satellite orbits. The columns give from left to right: the initial phase space coordinates of the satellite, 
x, y, z, v x , v y , v z ; the initial inclination to the Milky Way disc, i; the eccentricity, e; the sense of rotation; and the simulation output time. 
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Figure 4. The distribution of v z velocities at the solar neighbourhood, 8 < R < 9 kpc; \z\ < 0.35 kpc, for LMC-10°, LMC-20 and 
LLMC- 10°. The black line shows the Milky Way disc stars, the red the accreted satellite material. The distributions are normalised to 
peak at 1 and do not represent the mass in each component. Notice the prominent wings in the thin disc distribution for LMC-10° and 
LMC-20°: some of the stars originating in the thin disc show 'thick disc' kinematics. For the more massive LLMC-10°merger, the wings 
are now very prominent: the heated thin disc distribution looks like the accreted distribution in LMC-10 and little thin disc component 
remains. 



3.3 Results 

Our main results are shown in Figures [5] [7J and [9] Figure [5] 
shows the effect of increasing the impact angle; Figure!?] the 
effect of changing the satellite orbit; and Figure |5J the effect 
of changing the satellite and host galaxy properties. Some 
simulations were run for longer to ensure that the final state 
was in equilibrium. Typical times were ~ 5 Gyrs (see Table 

Figure |3] shows the distribution of v z velocities of stars 
at the solar neighbourhood for LMC- 10°, LMC-20 and 
LLMC- 10°. For LMC- 10° and LMC-20 , some of the stars 



originating in the Milky Way thin disc are scattered such 
that they have thick disc like kinematics (notice the promi- 
nent wings), but the majority of stars with thick disc kine- 
matics are accreted. For the more massive satellite merger 
- LLMC-10° - these wings are now very prominent and lit- 
tle thin disc remains. The heated thin disc distribution now 
looks like the accreted distribution in LMC- 10°. 

In most of what follows, we will decompose the final 
distributions into Milky Way stars and dark matter, and 
those accreted from the infalling satellite. This decomposi- 
tion makes sense, since we can expect the accreted satellite 
stars to be chemically distinct from the thin disc stars. Ob- 
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Figure 6. The satellite-disc inclination angle 9 as a function of 
time; lines are truncated when the satellite is fully accreted. For 
9 < 20° initially, the satellite is dragged into the disc plane; for 
9 = 40° initially, the satellite remains out of the disc plane at all 
times. Increasing the satellite mass (LLMC-10°) leads to faster 
disc plane dragging. 



servations at the solar neighbourhood show that thick and 
thin disc stars do show different chemistry, while there are 
so me thin disc stars w ith thick-disc like kinematics (3T]and 
see iBensbv et af1l2007ft . In 33X51 we investigate alternative 
decompositions where we include heated thin disc stars in 
the thick disc. 



3.3.1 Varying the impact angle 

The rows in Figure [S] show the effect of increasing the im- 
pact angle: simulations LMC-10°, LMC-20 , LMC-40 and 
LMC-60 . From left to right, the panels show: (a) logarith- 
mic density contours, viewed from side, in units of M q pc~ 2 ; 

(b) surface density as a function of J? in a slice \z\ < 1.1 kpc 
(the red dotted line for simulation LMC-10° shows an ex- 
ponential fit with scale length Ro = -R1/2/I.68 = 4.3 kpc); 

(c) density as function of z in a slice 8 < R < 9 kpc; (d) 
the rotation curve (v^(R) for \z\ < 0.35 kpc); and (e) the R 
(solid), cj> (dotted) and z (dashed) components of the stel- 
lar velocity dispersion as a function of projected radius. In 
all cases, we show the Milky Way stars (black), accreted 
satellite stars (red), Milky Way dark matter (blue dotted) 
and satellite accreted dark matter (blue). The black dashed 
lines (left three panels) and green lines (right panel) show 
the Milky Way disc initial conditions. The solid dots mark 
the half mass scale lengths for each component. 

For 9 < 20° we form a thick disc of stars and dark 
matter. At the solar neighbourhood (8 < R < 9 kpc) 
these lag the thin disc rotation (by 20-85 km/s); are hot- 
ter (o z ~ 40 — 60 km/s); and are of longer scale length and 



larger scale height {R1/2 ~ 7 kpc; Zi/ 2 ~ 0.9 kpc). The rota- 
tional lag of the stellar thick disc and a z both increase with 
the satellite inclination angle, 9. Notice that the thick disc 
stars lead the rotation towards the Galactic centre. This is 
because the thick disc distributions have a hole in the centre 
(see Figure Ob)). At the solar neighbourhood, a large frac- 
tion of thick disc stars are near apocentre and have small 
tangential velocity. Towards the Galactic centre, this cannot 
be maintained because of the decreasing density; most thick 
disc stars at radii ^ 4 kpc are near pericentre and have high 
tangential velocity. 

For 9 > 20° we no longer form a thick disc of stars or 
dark matter. Instead, we have boxy stellar and dark mat- 
ter distributions that are much hotter (a z > 100 km/s) and 
rotate more slowly (v c < 100 km/s). The boxy stellar mate- 
rial is a strong candida te for the new ' inner halo' structure, 
recently discovered by ICarollo et al. I l|2007h (see Table Q). 
A similar retrograde encounter could also explain some of 
the outer halo. Recall that these high inclination encounters 
(9 > 20°) are twice as likely as the lower inclination ones 
(©• 

The movement from disc-like to boxy structures is also 
clearly seen in Figure [5jc). Notice that for LMC-10°, there 
is a definite disc in the accreted stars (red line) and dark 
matter (blue line). For LMC-20 , the discs expand in scale 
height and length. For LMC-40 and LMC-6O , both distri- 
butions are constant in density over the range —3 < z < 3 
and can no longer be described as discs. 

Figure [6] demonstrates that dynamical friction plane 
dragging is responsible for the formation of accreted thick 
discs. In all simulations, the inclination angle of the satellite 
to the disc decreases with time. However, only for 9 < 20° 
initially, is the satellite completely dragged into the disc 
plane. 

3.3.2 Varying the orbit 

Figure[7]investigates the effect of changing the satellite orbit; 
lines and panels are as in Figure [S] We consider a retrograde 
encounter with the disc (LMC-Ret-10°); an orbit with low 
eccentricity (LMC-10°-eless; e = 0.36), and an orbit with 
large apocentre (LMC-10°-abig; r a — 81 kpc). In the cosmo- 
logical simulations of 321 we find no strong bias towards pro- 
grade or retrograde orbits, suggesting that retrograde merg- 
ers should be common (other recent stu dies also find onl y 
weak net rotation for the satellites, see e.g. lSales et alj2007h . 
From the observational point of view, it is interesting that, of 
the two extragalactic disc galaxies st udied in detail to date, 
one h as a counter-rotating thick disc (|Yoachim fc Dalcantonl 
2005). The orbital decay of the satellite for each of these runs 
and the reference simulation, LMC-10°, are shown in Figure 

m 

As compared to the equivalent prograde encounter 
(LMC-10°), the retrograde satellite disrupts further out 
leaving a larger hole behind. This can be seen in the ra- 
dial surface density (Figure Efb)). This is expected because 
dynamical friction goes as 1/v 2 . The prograde satellite is 
nearly stationary with re spect to the disc and falls in faster 
than the retrograde one (|Quinn fe Goodman1ll986l ). The ef- 
fect is small, however, because the dark halo dominates the 
mass as all radii in our MW model (this is not the case for 
MWB). The larger difference between the prograde and ret- 
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Figure 5. The effect of increasing impact angle: simulations: LMC-10°, LMC-20 , LMC-40 and LMC-60 . From left to right, the 
panels show: (a) logarithmic density contours, viewed from side, in units of Mq pc — 2 ; (b) surface density as a function of R in a slice 
\z\ < 1.1 kpc (the red dotted line [offset] for simulation LMC-10° shows an exponential fit with scale length Ro = iJj/2/1-68 = 4.3 kpc); 
(c) density as function of z in a slice 8 < R < 9 kpc; (d) the rotation curve (v<j,(R) for \z\ < 0.35 kpc); and (e) the R (solid), <f> (dotted) 
and z (dashed) components of the stellar velocity dispersion as a function of projected radius. In all cases, we show the Milky Way 
stars (black), accreted satellite stars (red), Milky Way dark matter (blue dotted) and satellite accreted dark matter (blue). The black 
dashed lines (left three panels) and green lines (right panel) show the Milky Way disc initial conditions. The solid dots mark the half 
mass scale lengths for each component. The black dotted lines mark the solar position. Notice that a thick disc of stars forms for impact 
angles ^ 20°, with a corresponding thick disc of accreted dark matter. Higher impact angles give more boxy stellar and dark matter 
distributions that rotate more slowly and are hotter. They also do increasing damage to the Milky Way thin disc, producing a flared 
outer disc that for LMC-40 and LMC-60 is better described as a warp. 



rograde encounters is in the response of the Milky Way disc; 
we discuss this in £13.3.41 

The lower eccentricity run, LMC-10°-eless (e = 0.36), 
gives a lower an and a smaller lag in the rotation of the 
accreted material. The small rotation lag can be increased 
by increasing the impact angle, 8 (see M3.3. II) . 

The high apocentre run, LMC-10°-abig, is particularly 
interesting. Even over 7Gyrs it has not fully accreted and 



dynamical friction proceeds slowly, despite the satellite be- 
ing initially of LMC mass. The orbit is shown in Figure [S] 
The dotted blue line on this plot shows the expected orbit 
for LMC-10°-abig assuming a static b ackground potential 
and C handrasekhar dynamical friction ( Binnev fc Tremaind 
2008). We use a Coulomb log arithm that varies with radius 
as in lHashimoto et al.l ([2003) . The simulated orbit decays 
much more slowly than the Chandrasekhar model because 
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Figure 7. The effect of changing the orbit: simulations: LMC-Ret-10°, LMC-10°-eless and LMC-10°-abig. The lines and panels are as in 
Figure[5] Notice that the retrograde merger, LMC-Ret-10°, forms a thick disc very similar to the equivalent prograde merger (LMC-10°; 
Figure^, though it disrupts further from the centre. It has a quite different effect on the Milky Way disc: a bar does not form, and there 
is a less pronounced warp/flare. The low eccentricity merger, LMC-10°-eless, produces a higher stellar density at the solar neighbourhood 
than our reference simulation LMC-10°, but the dark disc density is very similar. The high apocentre merger, LMC-10°-abig, fails to 
complete merging even over 7Gyrs (the remnant is visible in the left panel; see also Figure [8} and a significant thick disc is not formed. 



we have not accounted for satellite mass loss due to tides. 
The satellite, which starts out at 2.4 x 10 10 Mr, is just a few 
10 8 M at the end of the simulation. IZhaol (|2004T l have re- 
cently highlighted this problem. They show that it is difficult 
to understand how mergers - like the observed Sagittarius 
dwarf merger - could proceed in less than a Hubble time 
without Sagittarius having been unrealistically massive in 
the past. A solution, as presented in SJ21 is that galaxies like 
Sagittarius 'ride in' inside larger loosely bound groups and 
are effectively deposited at low apocentre. This is why in our 
other runs we assume an initial apocentre of r a — 30 kpc. 

3.3.3 Varying the satellite and host galaxy properties 

Figure [9] investigates the effect of varying the satellite and 
host galaxy properties. We consider a low mass satellite 
(Fornax-10°); a high mass satellite (LLMC-10°); and a more 
massive Milky Way disc (LLMC-10°-MWB). 

The Fornax merger (Fornax-10°) does not produce a 
structure resembling a thick disc. It disrupts too far out and 
is destroyed completely by R ~ 5 kpc. Its rotation velocity 
leads rather than lags the thin disc (because of the large hole 
interior to 5 kpc); its dispersions are too high; and the ac- 



creted stellar density is two orders of magnitude lower than 
LMC-10°. Summing over many such mergers could alleviate 
some of these discrepancies by filling in the hole with ac- 
creted material. However, nearly half of such mergers should 
be retrograde, while from we expect ~ 15 such mergers 
within 20° of the disc plane. The integrated light and dark 
matter from all such mergers will still be an order of magni- 
tude lower than that contributed by LMC-10 . We conclude 
that only the most massive mergers are of relevance for con- 
tributing accreted stars to the thick disc. 

While not forming thick discs, small-satellite encoun- 
ters do give rise to interesting fine-structure, in this case 
in the form of a ringjfl Since the satellite disruption radius 
will depend more on the satellite mass than on its orbit, 
ring features could originate from a super-position of satel- 
lites. Such super-positions will broaden the ring and make 



Rings can also be produced by low density satellites on more 
circular orbits, though such orbits seem unlikely in our current 
cosmology (see i[2l and Figure [TJ. High density satellites do not 
produce rings because they fall in via dynamical friction faster 
than they disrupt. 
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Figure 9. The effect of changing the satellite and host galaxy properties: simulations: Fornax-10°, LLMC-10° and LLMC-10°-MWB. 
The lines and panels are as in Figure[5] The low mass satellite, Fornax-10°, does not perturb the Milky Way disc at all, nor does it form 
anything which may be reasonably described as a thick disc. The high mass satellite, LLMC-10° produces results similar to LMC-10°, 
but docs significant damage to the Milky Way disc which now has density comparable to the dark halo at the solar neighbourhood. 
The more massive Milky Way disc run, LLMC-10°-MWB, produces results similar to LLMC-10°, but the increased disc mass causes a 
prominent bar to form, while the satellite disrupts more rapidly leading to a thick disc of larger scale height and lower density interior 
to 2 < 1.1 kpc. 



it appear disc-like at large radii. Indeed, this could provide 
an explanation fo r the extended disc seen in Andromeda 
|lbata et al.ll2005l and flfl From [H we expect ~ 15 such 
low inclination, low mass, mergers to occur in the Milky 
Way. Some of these will have occurred before the disc of the 
Milky Way was fully formed and may be essentially unde- 
tecable today. Howe ver, one such stru cture does appear to 
have been observed (|Conn et al.ll2007h . There may be more 
to be found. 

LLMC-10° investigates the effect of increasing the satel- 
lite mass. Many of the properties of the accreted mate- 
rial are reminiscent of our reference run - LMC-10°: it 
lags the rotation of the thin disc by ~ 40km/s, is hot- 
ter with a z ~ 80km/s at the solar neighbourhood and 
is of larger scale height (zi/2 = 0.9 kpc) and longer scale 
length (Ri/2 = 8 kpc) than the thin disc. Such mergers - 



9 Alternative explanations in the literature are to have on e large 
satellite merge in the disc-plane l lPenarrubia et alJlioO^) : or to 
have a prograde fly-by from a massive satellite o n an eccentric 
orbit llKazantzidis et al.ll2007l ; lYounger et al.ll2008l) . 



Wmax — 80km/s - are common; we expect on average one 
such merger to occur within 20° of the disc plane (assum- 
ing isotropy). Furthermore, what is particularly interesting 
about this simulation is the heating of the thin disc. The 
final thin disc scale height and kinematics for this run - un- 
like all of the previous runs - now matches the thick disc of 
the Milky Way with Z\ii = 0.6 kpc and a z — 40 kpc at the 
solar position. We will discuss this in more detail in i|3.3.4l 
and 33331 



Finally, LLMC-10°-MWB investigates the effect of in- 
creasing the mass of the Milky Way disc. With the more 
massive Milky Way disc, the satellite now disrupts further 
out than in LLMC-10 , leading to a thick disc of larger scale 
height and lower density interior to \z\ < 1.1 kpc. The more 
massive disc goes bar unstable without any encounter, form- 
ing a prominent bar that appears bulge-like when viewed 
from the side in project ion (a detailed account of this type 
of instability is given in iDebattista et aljfeood ). 
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Figure 8. Decay of the satellite orbit as a function of time for 
selected simulations. The dotted blue line shows the analytic de- 
cay rate expected for LMC-10°-abig, assuming Chandrasekhar 
friction and that the satellite does not lose any mass. 



3.3.4 The Milky Way thin disc: a bar, flare, warp and 
heated thick disc 

Figures 1101 1111 and the left-most panels of Figures [5] [7] and 
show the effect of the mergers on the underlying Milky 
Way disc. These results should be treated with some caution 
since we do not know what the Milky Way disc looked like 
at z = 1 and we have assumed near-present day properties 
for our initial conditions. Nonetheless, the mergers produce 
a wealth of interesting substructures in the disc that match 
observations of our own Galaxy, and are long-lived (recall 
that these simulations have been evolved for ~ 5 Gyrs or 
longer) . 

Figure [TU] shows the half mass scale height of the Milky 
Way disc z x /i as a function of R. The different coloured 
lines show sel ected simulations; the d a shed line shows a fit 
to data from lLopez-Corredoira et alj ((2002) for the Milky 
Way flare: 

Z1/2 (R) = 3.6 x ICT^e e X p ( 12 _ fl 6fl g Q pc)kpc ) (4) 

with i? = 7.9 kpc. 

The first peak in the distribution is due to the bar, 
the secondary rise is due to the flared disc. Increasing the 
impact angle makes very little difference to the size of the 
bar in scale length or height and in all cases we obtain 
Z\/2 ~ 0.3 kpc; R ~ 2 kpc. (Changing the satellite orbit also 
makes almost no discernible difference and we omit these re- 
sults for clarity.) Increasing the satellite mass increases the 
bar size by only a modest amount (R ~ 2.5 kpc). This is 
interestin g because the Milky Way bar a ppears to be longer 
than this. lLopez-Corredoira et alj (|2007T ) have recently con- 
firmed the existence of the Milky Way bar using 2MASS star 



counts. They find dimensio ns of 3.9 x 1.2 x 0.2 kpc. Similarly, 
iBissantz fc Gerhard! (|2002h have recently found a bar length 
of 3.5 kpc. The only bar we form that is this long is in LLMC- 
10°-MWB (R ~ 4 kpc). This bar forms as a result of secular 
evolution in the disc and is not driven by the merger. It is 
also significantly fatter (21/2 = 0.9 kpc) than the Milky Way 
bar, though many of these stars would likely be classified as 
bulge, rather than bar starfj. A final interesting point is 
that the retrograde orbit LMC-Ret-10° hardly excites a bar 
at all. S uch results have been seen in the literature dating 
back to iHolmb erd (|l94lf ). who presente d the first study of 
prograde versus retrograde interactions. iToomre fc Toomrel 
|l972T ) attribute the phenomenon to increased resonance in 
the prograde case, but one can also think of it in terms of 
increased tidal for ces for the prograde interaction (see e.g. 
iRead et alj|2006bl ). 

The flare and warp are much stronger indicators of 
merger activity than the bar. From Figure 1101 we can see 
that the strength of the flare grows with satellite impact 
angle and is not affected by any other parameter - particu- 
larly at large radii (R > 10 kpc). LMC-6O provides the best 
match to the observed Milky Way flare. This is encourag- 
ing since such mergers should be twice as likely as the low 
inclination mergers that produce the accreted thick disc. 

In Figure [11] we show the mean height of the disc along 
x' - the direction of maxim al height variation; the dashed 
line shows a fit to data from lLopez-Corredoira et al.1 (|2002l ) 
for the Milky Way warp: 

< z >= 1.2 x 10 _6 i?(kpc) 5 ' 25 sin(<£ + 5°)kpc (5) 

We average over \y'\ < 0.5 kpc and all z. As with the flare, 
only LMC-6O provides a good match to the observed Milky 
Wa y warp. 

iLevine et al.1 {2006) have recently suggested that the 
Milky Way warp is a result of resonant driving by the LMC, 
but it is intriguing that high impact mergers can produce 
similar effects, while simultaneously providing a good match 
to the Milky Way flare. A natural candidate for producing 
the warp is t he Sagittarius dwarf which recently fe ll in on 
a polar orbit IjBailin 2003; B ailin fc Steinmetal2005l ). How- 
ever, since such merger produced warps are long-lived, ear- 
lier mergers could have been more important. 

The longevity of th e warp is perhaps surprising. 
iNelson fc Tremainei (|l995h calculated that dynamical fric- 
tion ought to rap idly dam p such warps away (se e also 
iBinnev et aHll998f >. However. IShen fc Sellwoodl (2006) show 
that this does not occur in practice because the in- 
ner halo moves with the disc. They found, like us, long 
lived warps that survive ove r several Gy rs (see also 
iDebattista fc Sellwoodlll999l and lBinney||2007l for a review). 

Finally, from Figure[TU]notice that the two largest merg- 
ers (LLMC-10°and LLMC-10°-MWB; tw ~ 80km/s) heat 
the thin disc enough that it resembles the thick disc of the 
Milky Way (21/2 (ii©) ~ 0.6 kpc). None of the lower mass 
satellites produce enough vertical heating. We discuss the 

10 The half mass scale height for the Milky Way bulge quoted 
in Table [T] is 0.25 kpc which seems much smaller than the 0.9 kpc 
for the secular bugle formed in LLMC-10°-MWB. However, the 
difference is not so acute. The data for the Milky Way bulge 
probe only the central ~ 1 kpc, where our secular bulge has a 
comparable Z1/2 ~ 0.3 kpc. 
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Figure 10. The half mass scale height of the Milky Way disc 
2]/2 as a function of R. The different coloured lines show se- 
lected simulations; the dashe d line shows a fit to data from 
lLopez-Corredoira et all [120021) for the Milky Way flare. The first 
peak in the simulated distributions is due to the bar, the sec- 
ondary rise is due to the flared disc. 
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Figure 11. The mean height of the disc along x 1 - the direction of 
maximal height variation. We average over \y'\ < 0.5kpc and all z. 
The different coloured lines sho w selected simulations; t he da shed 
line shows a fit to data from lLopez-Corredoira et al. (20q2) for 
the Milky Way warp. 



contribution of these heated thin disc stars to the thick disc 
in 33~331 

Our results for the response of the Milky Way disc 
agree very well with othe r recent studies in the literature. 
iDubinski et"ail (|2008h and lGauthier et al.1 (|2006h merge 100 
satellites with a simulated M31 galaxy. They find, like us, 
a warp, bar and flare that are ex cited by LMC-mass inter- 
actions. iKazantzidis et all (|2007T > also find a warp, bar and 
flare form naturally as a result of disc-satellite interactions 
in a ACDM cos mology. Both of these stud ies are comple- 
menta ry to ours. IDubinski et alj |2008h and lGauthier et al.l 
(2006) consid er the effect of many m ergers that occur si- 
multaneously, IKazantzidis et all ((2007) consider the cumu- 
lative effect of many interactions that each pass the disc 
only once, and we follow single interactions over many or- 
bits until they merg e. An important dif f erence between 
IDubinski et~ail l|200Sli and iGauthier et~ail ^OOrj ) and our 
study is that they merge the present surviving satellite distri- 
bution with M31 (this is al so the case for mos t previ ous stud- 
ies, e.g. iFont et~aLll200ll '). IKazantzidis et al.1 (|2007m . like us, 
use the distribution of accreted satellites that is systemati- 
cally more massive and more destructive. As a result, these 
earlier studies underestimate the vert ical heating of the disc . 
A se cond important differe nce is that IDubinski et al.l (|2008l ) 
and IGauthier et~aH (|2006T ) do not account for the larger 
group environment around the subhalos. They find that dy- 
namical friction is largely unimportant, as in our LMC-10 - 
abig simulation. Once the larger group environment is taken 
into account, dynamical friction does play an important role 
and allows massive subhalos to successfully merge with the 
disc at high redshift. Finally our results for retrograde merg - 
ers agree well with a study by IVelazquez fc White! l|l999l) . 
They found, like us, that retrograde interactions give very 
little disc heating. This demonstrates that the heating is a 
result of resonant driving from the satellite that is strongly 
suppressed in retrograde interactions. 

3.3.5 Accreted versus heated thick discs 

Combining the results from the previous sections, we find 
that an accreted thick disc of stars forms for impact angles 
^ 20°. For an impact angle of 10° < 8 < 20°, an eccentric- 
ity of 0.36 < e < 0.8, and a satellite mass of 70 — 90km/s, 
the accreted material has almost all of the observational 
properties of the Galactic thick disc at the solar neighbour- 
hood (8 < R < 9kpc). For our reference simulation, LMC- 
10°, we find a thick disc that lags the thin disc rotation 
by ~ 20km/s, is hotter with a z ~ 40km/s, compared with 
a z ~ 20km/s for the thin disc, and has a larger scale height 
(zx/2 = 0.7 kpc) and longer scale length (R1/2 — 7.3 kpc; 
Ro — 4.3 kpc) than the thin disc. Increasing 9 increases a z 
and the rotation lag; reducing e reduces or and the rotation 
lag; and increasing the satellite mass increases the velocity 

dispersion of the thick disc in all directions. 

However, as in previous studies (e.g. [Walker et al.l 

1996), our accreted thick discs are significantly less mas- 
sive than that of the Milky Way. In LMC-10°, we find 
Pthick/pthin = 0.4% (assuming the observed value of 
Pthin = 0.09 Mq pc" 3 ), which is too low by a factor 
~ 30 (see Table [TJ . Our other runs produce similar num- 
bers: LLMC-10° gives the most massive thick disc, with 
Pthick/pthin = 1.3%, but even this is too small by a factor 
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Figure 12. Current constraints on accreted stellar and dark 
discs at the solar neighbourhood. The plot shows the dark disc 
to dark halo ratio Pddisc/phalo as a function of the thick to 
thin disc density ratio Pthick/pthin; we assume for this plot 
PHALO = 0.01 M pc -3 and pthin = 0.09 M pc -3 . Results 
are extrapolated from two simulations: LMC-10° and LLMC- 
10°, marked by the diamond and square, respectively. The star 
and plus show two additional simulations: LLMC-10°-MWB and 
LLMCB-10°-MWB, respectively. The different colours extrapo- 
late to lower mass to stellar mass ratio, M/M„. The solid lines ex- 
trapolate to multiple mergers of identical satellites. The X's show 
the effect of including heated thin disc stars in the thick disc for 
LMC-10° (60%) and LLMC-10° (90%; see text for details). The 
vertical dotted black line marks the observed Pthick/pthin at 
the solar neighbourhood. 



~ 10. In fact, these normalisations are a much close r match 
to the local phalo/pthin = 0.5% (jjuric et al.ll2008l ). 

The thick disc density normalisation Pthick/pthin can 
be increased by: (i) increasing the satellite stellar mass; (ii) 
integrating over many near disc plane mergers; and (iii) in- 
cluding heated thin disc stars in the thick disc. A possible 
fourth mechanism for driving up the mass at the solar po- 
sition would be to put the stars and dark matter in a ring, 
rather than a disc. In practice this does not work because dy- 
namical friction drives the satellite inwards as fast as it loses 
mass and the resulting density at the solar neighbourhood 
remains low. (We tested this explicitly using test simulations 
where the satellite was placed on a circular orbit; these are 
omitted for brevity.) 

In Figure 1121 we summarise each of the effects (i) 
to (iii), above by extrapolating the results for LMC-10 
and LLMC-10°. Our extrapolations were performed in the 
following way: 

(i) Increasing the satellite stellar mass 

We increase the satellite stellar mass while holding the 
dark matter mass fixed. This is shown by the black red 



Figure 13. The mass in stars (red) and dark matter (blue) that 
contribute to the thick discs interior to \z\ < 1.1 kpc, as a func- 
tion of their initial radius within the satellite (in 1 kpc bins) for 
LMC-10° (solid lines); LLMC-10° (dotted lines); and LLMCB- 
10°-MWB (dashed lines). 



and blue lines in Figure [121 which are for M/M r = 35, 15, 6, 
respectively. 

To test the extrapolation, we ran an additional simula- 
tion LLMCB-10°-MWB (shown by the star in Figure US). 
This was identical to LLMC-10°-MWB (shown by the plus) 
except that the satellite had twice the stellar mass. The 
mass in accreted stars doubled, as expected. However, the 
mass in the dark disc also increased. The reason for this 
is given in Figure 1131 This shows the mass in stars (red) 
and dark matter (blue) that contribute to the thick discs 
interior to \z\ < 1.1 kpc, as a function of their initial radius 
within the satellite for LMC-10 (solid lines); LLMC-10 
(dotted lines); and LLMCB-10°-MWB (dashed lines). Most 
of the mass that contributes originates from the central 
10 kpc in the satellite. This is the majority of the stars, but 
only a fraction of the dark matter. For LLMCB-10°-MWB, 
the increased disc mass caused the central satellite dark 
matter to adiabatically contract, causing more dark matter 
to end up in the thick disc too. This test demonstrates that 
it is the satellite mass interior to ~ 10 kpc - both stars and 
dark matter - that contributes to the thick disc. It also 
explains why the mass to light ratio of our accreted thick 
discs is lower than that of the satellites from which they 
form. 

(ii) Multiple mergers of like satellites 

The solid lines in Figure Q2] show the effect of multiple 
mergers of like satellites. We assume that the accreted 
material simply adds to that already present. From Figure 
[1] left panel, we can see that cosmic variance gives an 
upper bound of ~ 10 mergers above LMC mass. However, 
for multiple mergers to contribute to the thick disc, they 
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must all be low inclination and prograde. Such a scenario 
seems unlikely unless accretion occurs along filaments and 
the disc is favourably aligned. 

(Hi) The effect of heated thin disc stars 
For our most massive mergers - LLMC-f 0° and LLMC-f 0°- 
MWB - the originally thin disc is heated sufficiently that it 
resembles the Milky Way thick disc at the solar neighbour- 
hood: or ~ 80km/s; a$ ~ 60km/s; a z ~ 40km/s; rotation 
lag ~ 20km/s; and scale height zy 2 ~ 0.6kpc. 

Here we consider the effect of adding these heated thin 
disc stars to the thick disc, both for the most massive merg- 
ers, where almost all of the stars can contribute, and for the 
less massive merger LMC-fO°. We populate the thick disc 
with thin disc stars at the solar neighbourhood by compar- 
ing thick and thin disc distribution functions for each star: 

v 2 ("^-<^>) 2 „2 2 

f = Ae~*°?e e~^?e~^ (6) 

with o r — 85km/s, o<$, = 60km/s, o z — 40km/s, < v$ >= 
f83km/s, and z e = f.48kpc for the thick disc; and o r = 
60km/s, (70 = 40km/s, a z — f9km/s, < v# >= 200km/s, 
and z e = 0.85 kpc for the thin disc. These choices are mo- 
tivated by our results in H3.3.1I but our results are not sen- 
sitive to these parameters. We assign thin disc stars to the 
thick disc if /thick//thin > 1 for that star. We then vary 
the fraction of thin disc contaminants by varying the relative 
normalisation of each distribution, A. 

For our reference simulation, LMC-fO°, we find that 
we can populate up to 60% of the thick disc with hot thin 
disc stars without significantly reducing the thick disc scale 
height. However, greater contamination than this leads to 
a thick disc that is too cold vertically. For LLMC-f 0°, the 
more massive satellite heats the thin disc enough that it re- 
sembles a thick disc. A decomposition is still required, how- 
ever, since the heated thin disc at the end of the simula- 
tion does not show enough rotation lag (see Figure [9jd)). 
We find that for this simulation 90% contamination gives a 
good qualitative match to the Milky Way thick disc and re- 
covers the full observed thick disc mass. We show the effect 
of 60% thin disc contamination for LMC-f0° and 90% for 
LLMC-f 0° by the X's in Figure [H 

From Figure If 21 we can see that while multiple merg- 
ers and lower M/M* help to increase the thick disc mass, 
heated thin disc stars due to a very massive merger (v max ~ 
80 km/s) seem essential for reaching the high thick disc mass 
observed in the Milky Way. Accreted stars likely contribute 
fO — 50% of the observed thick disc. This is in reason- 
able agreem e nt wit h the recent cosmological simulation of 
lAbadi et al.l (|2003h . where they find a thick disc that is 
equally composed of heated thin disc stars and accreted 
material. IVillalobos fc Helmil (|2008l ). also find that accreted 
stars can make up only a small fraction of the thick disc in 
the disc plane. 

The above seems at odds w ith the recent extra- 
galac tic observations of FGCf4f5 (jYoachim fc Dalcantonl 
120051 ). which shows a massive counter-rotating thick disc. A 
possible solution is that FGC14f5 had a rather more mas- 
sive (and more rare) retrograde merger than those studied 
here. 

In principle, heated and accreted stars should be sepa- 
rable by their chemistry. The high Pthick/pthin obtained 



at the solar neighbourhood comes from a fit to star counts 
for stars that are well separated in colour, bu t for which 
there is currently no spectroscopic information l|juric et al.l 
120081 ). But this work will be difficult. State-of-the-art spec- 
troscopic studies have too fe w stars (~ 200) to obtain an 
accurate density distribution (|Bensbv et al.ll2007l ). 

As a final note for this section, recall that the Milky 
Way thick disc has stellar mass ~ I0 10 Mq, while the stel- 
lar halo is ~ 10 9 Mq (see Table [IJ. If we put too many 
accreted stars into the thick disc then we have a problem: 
the stellar halo will become too massive because there will 
be two high inclination mergers for every merger near the 
disc plane. This suggests that either the Milky Way had its 
most massive merger near the disc plane, or - which is more 
likely - a significant fraction of the Milky Way thick disc 
comprises heated thin disc stars (or the thick disc formed 
through some other mechanism). 

3.3.6 A dark matter thick disc 

A key new idea presented in this work is that near-disc ac- 
cretion - that must occur in a ACDM cosmology - leads 
to the formation of a dark matter disc, as well as an ac- 
creted stellar thick disc. We can see this in the left three 
panels of Figure [5] (blue lines). For impact angles < 20°, a 
dark disc forms. For LMC-f 0° it has scale length and height 
given by: R1/2 = ff.7kpc, z\/2 = 1.5 kpc, while its den- 
sity at the solar neighbourhood is Pddisc (^0) = 0.22phalo 
(assuming Phalo = 0.01 M pc -3 ). For LLMC-10° it has 
Z1/2 = 1.7kpc and pddisc(-Rg) = 0.42p H ALO- At these low 
densities, the dark disc is not dynamically interesting in 
the disc plane. Integrating over the total expected merg- 
ers brings the density up, but it is hard for it to exceed the 
density of the dark halo at the solar neighbourhood without 
having an unrealistically high number of mergers (see JZJ. 
The extrapolated constraints for LMC-10° and LLMC-10° 
are shown in Figure If 21 

We can obtain an upper bound on pddisc for \z\ < 
1.1 kp c from the kinematics of s tars at the solar neighbour- 
hood. |Hgtaberg^&j^]xnn| (120001 ) find a local dynamical mass 
density of 0.102 ± 0.01M Q pc" 3 with 0.095 M pc -3 in visi- 
ble disc matter. This leaves room for a maximal dark density 
of 0.008 M pc -3 which seems consistent with a da rk halo 
and n o dynamically significant dark disc. However, IStatlerl 
( 1989) have argued that model degeneracies in this analysis 
could lead to systematic errors as large as 30%. If we take 
this more generous limit, we can have a maximal density of 
0.038 Mq pc -3 , which leaves room for a dynamically signif- 
icant dark disc up to Pddisc (Rq) = 3phalo- These limits 
will improve dramatically with the improved data from large 

kinematic surveys l ike RAVE and G AIA. 

Interestingly, IfCalberlal (|2003h and iKalberla et al.l 

(2007) have recently argued for a dark disc in the Milky Way 
to explain the observed flaring in the HI gas. The dark disc 
they require has a scale height and length of Z1/2 ~ 2.8 kpc, 
R 1/2 = 12.6 kpc and total mass ~2x 10 11 Mq. For LLMC- 
10°, our dark disc has total mass 8.6 x 10 9 for \z\ < 1.1 kpc 
and 3 x 10 10 Mq for \z\ < 5 kpc; theirs has 4 — 6 times more 
mass than this within the same z-limits. While the half mass 
sc ale length and height of our dark disc match well those 
of IKalberla et all l|2007h , our disc appears significantly less 
massive. 
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Figure 14. Distribution functions for the thin disc (black), thick disc (red), dark disc (green) and dark halo (blue) for LMC-10°at the 
solar neighbourhood (8 < R < 9kpc; \z\ < 0.35 kpc); we plot results from the Earth frame so that the dark halo appears to rotate, while 
the discs are stationary. The distributions are normalised to peak at 1 and do not represent the mass in each component. Notice that 
the dark disc has very similar kinematics to the stellar thick disc, except in that it has slightly larger a z . Similar results were found for 
our other runs. The radial velocity distribution is triple-peaked because the stars and dark matter are not fully phase mixed yet, even 
after 6.5 Gyrs. 



Even if dynamically uninteresting, the dark disc has im- 
portant implications for the direct detection of dark matter 
because it rotates much more slowly with respect to the 
Earth than the Milky Way halo. The solar neighbourhood 
distribution functions for the thin disc, thick disc, dark disc 
and dark halo in LMC-10° are shown in Figure ITU we plot 
results from the Earth frame so that the dark halo appears 
to rotate, while the discs are stationary. Notice that the dark 
disc has very similar kinematics to the stellar thick disc. If 
the thick disc of the Milky Way is mostly accreted rather 
than heated stars, then this is useful. It means that we can 
obtain a good measure of the kinematics of the dark disc 
from observations of the stellar thick disc at the solar neigh- 
bourhood. 

The effect of the dark disc on the dark matter detection 
rate, the annual modulation signal, and on the dark matter 
capture rate by the Sun an d Earth will be qu antified in de- 
tail in forthcoming papers (jBruch et al . 2008; Bruch, Read, 
Baudis and Lake in prep. 2008). 



4 CONCLUSIONS 

We have presented a detailed study of the effect of the Milky 
Way stellar disc on merging satellites in a ACDM cosmol- 
ogy. We have used both cosmological dark matter simula- 
tions of structure formation to assess the likelihood of near- 
disc plane accretion events, and collisionless simulations of 
satellite mergers to quantify the final state of the accreted 
material, and the effect on the Milky Way stellar disc. Our 
main findings are as follows: 

(i) On average, a Milky Way-sized galaxy has 3 sub- 
halos with peak circular speed « max > 80km/s; 7 with 
«max > 60km/s; and 15 with f max > 40km/s merge at 
redshift z 1. Assuming isotropy, a third of these merge at 
an impact angle 8 < 20° and are dragged into the disc plane 
by dynamical friction. 

(ii) For an impact angle of 10° < < 20°, an eccentricity 



of 0.36 < e < 0.8, and w max = 60 — 80km/s, the accreted 
material has almost all of the observational properties of 
the Galactic thick disc at the solar neighbourhood (8 < R < 
9 kpc), but the resultant accreted thick disc is under-massive 
by a factor ~ 2 — 10. 

(iii) The most massive mergers at high redshift - 
u max ^ 80km/s - heat the thin disc enough to produce a 
heated thick disc; none of the lower mass mergers provide 
enough heating. The heated thick disc is essential for obtain- 
ing a thick disc as massive as that observed in the Milky Way 
and likely contributes ~ 50 — 90% of the thick disc stars at 
the solar neighbourhood. 

(iv) A key new point we make is that low inclination 
mergers - that must occur in a ACDM cosmology - also give 
rise to a thick disc of dark matter. The dark disc is of longer 
scale length and height than the accreted stellar thick disc 
(Ri/2 ~ 12 kpc, z\/2 ~ 1.5 kpc), and provides ~ 0.25 — 1 
times the density of the dark halo within |z| < 1.1 kpc at 
the solar neighbourhood. The precise number depends on 
the satellite properties and the number of mergers. At all 
but the highest of these densities, the dark disc is not likely 
to be dynamically interesting. However, the dark disc does 
have important implications for the direct detection of dark 
matter because of its low velocity with respect to the Earth. 
A full quantitative ana lysis is presented in a companion pa- 
per jBruch et al.ll2008r i. 

(v) Higher inclination encounters 6 > 20° are twice as 
likely as low inclination ones. These lead to structures that 
are hotter [a z > lOOkm/s), rotate more slowly (v c < 
lOOkm/s), extend 10 — 20 kpc above the disc, and closely 
res emble the inner /ou ter stellar halos recently discovered 
bv lCarollo et aL (2007). The integrated light from such en- 
counters is of order the total mass in the Milky Way stellar 
halo. 

(vi) We quantify the effect these mergers have on the 
Milky Way disc. All encounters excite a bar, flare and warp 
in the disc that are long-lived. The bar is a poor indicator of 
merger activity because a bar formed by secular evolution 
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(due to a more massive stellar disc) is longer and provides a 
better fit to the Milky Way bar than those induced by merg- 
ers. By contrast, the flare and warp are strong indicators of 
merger activity. Their strength grows with satellite impact 
angle and is not strongly affected by any other parameter - 
particularly at large radii (R > lOkpc). Of the simulations 
that we consider, = 60km/s merger at 9 = 60° to 

the Milky Way disc plane provides the best match to the 
observed Milky Way warp and flare. 

(vii) Finally, we show that even w m ax = 70km/s galax- 
ies have dynamical friction times longer than the age of 
the Universe once mass loss due to tides is taken into ac- 
count. In practice, in our cosmological simulations, they 
merge much faster than this because they 'ride in' inside 
loosely bound groups that are ten times more massive. This 
should be taken into account when calculating the orbits of 
Local Group satellite galaxies such as the Large and Small 
Magellanic Clouds, and the Sagittarius dwarf galaxy. 
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